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Abstract 

Purpose Halting the loss of biodiversity while providing food 
security for a growing and prospering world population is a 
challenge. One possible solution to this dilemma is organic 
agriculture, which is expected to enhance biodiversity on the 
farmland. However, organic products often require larger 
areas. This study demonstrates how we can quantify and 
compare the direct land use impacts on biodiversity of organic 
and conventional food products such as milk. 

Material and methods This study assessed direct land use 
impacts of 1 1 of milk leaving the farm gate. Inventory data 
on land occupation were extracted from a life cycle assess¬ 
ment study of 15 farms in southern Sweden. Direct land use 
change data were derived from the FAO statistical database. 
Spatially differentiated characterization factors of occupation 
(CF 0cc ) and transformation (CF Trans ) were calculated based on 
the relative difference of plant species richness on agricultural 
land compared to a (semi) natural regional reference. Data on 
plant species richness and regeneration times of ecosystems 
(for calculating transformation impacts) were derived from a 
literature review. To account for differences in biodiversity 
value between regions, a weighting system based on absolute 
species richness, vulnerability and irreplaceability was 
applied. 
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Results and discussion Organic milk had a lower direct land 
use impact than conventional milk, although it required about 
double the area. Occupation impacts dominated the results and 
were much smaller for organic than conventional milk, as 
CFqcc of organic land uses were considerably smaller. For 
transformation impacts, differences between the two farming 
practices were even more pronounced. The highest impacts 
were caused by soymeal in concentrate feeds (conventional 
milk) due to large-scale deforestation in its country of culti¬ 
vation (i.e. Brazil and Argentina). However, lack of reliable 
data posed a challenge in the assessment of transformation 
impacts. Overall, results were highly sensitive to differences 
in land occupation area between farms, the CF 0cc and assump¬ 
tions concerning transformed area. Sensitivity and robustness 
of results were tested and are discussed. 

Conclusions Although organic milk required about twice as 
much land as conventional, it still had lower direct land use 
impacts on biodiversity. This highlights the importance of 
assessing land use impacts not only based on area but also 
considering the actual impacts on biodiversity. The presented 
approach allows to quantify and compare hot- and coldspots 
in the agricultural stage of milk production and could poten¬ 
tially also be applied to other agricultural products. However, 
more research is needed to allow quantification of indirect 
land use impacts. 

Keywords Biodiversity • Land use • Life cycle impact 
assessment • Milk • Organic agriculture 

1 Introduction 

One of the most important drivers of global biodiversity loss is 
land use change or habitat destruction (Sala et al. 2000; Foley 
et al. 2005; Pereira et al. 2010), which is indirectly linked to 
increased food and energy demand by a growing and prosper¬ 
ing world population (Tilman 2001). Providing food security 
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for the expected future human population of 9 billion people 
(Zhang 2008) while halting the loss of biodiversity is a chal¬ 
lenge. Over the past years, there has been a growing debate 
about whether land sparing (i.e. agricultural intensification to 
minimize pressure to convert wild nature) or wildlife-friendly 
farming (i.e. increasing biodiversity on farmland which might 
reduce agricultural yields) is the more promising solution 
(Tschamtke et al. 2012; Chappell and LaValle 2011; Ewers 
et al. 2009). Organic farming practices are found to increase 
species abundance (Crowder et al. 2010; Rahmann 2011) and 
species richness of many taxa (Bengtsson et al. 2005; Fuller 
et al. 2005; Gomiero et al. 2011) as well as the presence of 
endangered species (van Elsen 2000) compared to conven¬ 
tional farming. This positive effect on biodiversity points to 
the management principles of organic agriculture, which pro¬ 
hibit the use of chemical fertilizer or inorganic pesticides 
(Hole et al. 2005) and are characterized by a lower stocking 
rate (Hansen et al. 2001) and a more varied crop rotation 
(Maeder et al. 2002). 

Especially in Europe, efforts to preserve the high biodiver¬ 
sity value of long-established agricultural ecosystems are high 
(Green 2005). Therefore, subsidies for organic farming were 
introduced in the early 1990s (Stolze and Lampkin 2009). In 
2005, the European Union’s (EU plus national funds) support 
of organic farming increased to 17 % of the yearly expendi¬ 
tures (or a total of 0.6 billion Euro) dedicated to agri¬ 
environment measures (European Union 2007). In order to 
assess the usefulness of these large amounts of subsidies for 
biodiversity conservation, the overall impacts of organic and 
conventional products on biodiversity need to be quantified. 

An important decision support tool in this context is life 
cycle assessment (LCA), which allows the quantification of a 
range of environmental impacts throughout all stages of a 
product’s (often globally distributed) life cycle (e.g. Mila i 
Canals et al. 2007). Several approaches towards the assess¬ 
ment of direct land use impacts on biodiversity have been 
published. However, these studies are either only applicable to 
one or two world regions, and thus cannot capture land use 
impacts occurring in many world regions (e.g. Koellner 2000; 
Koellner and Scholz 2007; Michelsen 2008; Schmidt 2008a), 
or do not capture differences among different farming systems 
(e.g. Weidema and Lindeijer 2001; de Baan et al. 2013; 
Mattila et al. 2012). 

In this study, we quantify and compare the direct land use 
impacts on biodiversity of organic and conventional milk. We 
therefore further develop the methodology presented in de 
Baan et al. (2013) and apply it to a case study of milk 
production in Sweden, where organic milk was shown to 
require considerably larger areas compared to conventional 
milk (Cederberg and Flysjoe 2004). Milk is an interesting 
example because its land use is distributed over many world 
regions via the supply chain of concentrate feed (Cederberg 
and Mattson 2000). Thus, a regionalized global approach (as 


suggested in de Baan et al. 2013) is needed to accurately 
assess biodiversity impacts, due to its heterogenous distribu¬ 
tion and variable response to land use. In addition, the com¬ 
position and origin of feedstocks for organic and conventional 
milk vary substantially. For producing organic milk in the EU, 
more than 60 % of the daily intake of cows has to be roughage 
such as silage, hay or pasture and more than 50 % of the feed 
has to be produced on the farm itself (EC 2007). Conventional 
dairy cows are fed more on concentrate feed compared to 
organic cows, and the feedstock are to a larger degree 
imported from sub-/tropical countries (Cederberg and Mattson 
2000). In these countries, large areas of natural vegetation of 
high species richness and endemism (Brooks et al. 2006; Kier 
et al. 2005; Myers et al. 2000) have been transformed in the 
last decades (FAOSTAT 2012a; Mayaux et al. 2005). 

Milk is not only an interesting case to illustrate the trade¬ 
offs between area requirement and impacts on biodiversity but 
also an important commodity. The consumption of milk prod¬ 
ucts is projected to rise considerably in the next decades 
(Delgado 2003). Many LCA studies have been conducted on 
the comparison on organic and conventional milk. In terms of 
energy use, eutrophication and global warming potential no 
clear result of organic compared to conventional milk can be 
found (Thomassen et al. 2008; Flysjoe et al. 2012). In contrast, 
land use (in terms of area) was found to be considerably larger 
for organic milk than for conventional milk (Cederberg and 
Mattson 2000; Thomassen et al. 2008; van der Werf et al. 
2009). However, the impacts of land use on biodiversity have 
not yet been incorporated into LCA studies of milk (Yan et al. 
2011 ). 

2 Material and methods 

2.1 Description of compared agricultural systems 

Data on land occupation area and country of cultivation of 
fodder crops were taken from a study on dairy farms in south¬ 
west Sweden (Cederberg and Flysjoe 2004), following the 
same approach as Cederberg and Mattson (2000). In this 
study, we considered nine high-intensity conventional farms 
and six organic farms. In terms of total farm area, the two 
agricultural systems were comparable; however, the conven¬ 
tional farms had to purchase more concentrated feed as they 
had more cows on their farms (Table 1). The higher share of 
roughage feed on the diets of organic cows results in lower 
milk yields of organic compared to conventional cows. The 
functional unit (FU) of the study was " 1 kg of energy corrected 
milk" leaving the farm gate, i.e. transportation and processing 
of raw milk were excluded (Cederberg and Flysjoe 2004). 
Partitioning of environmental impacts (i.e. allocation) be¬ 
tween meat and milk as well as co-products in purchased 
concentrate feed, such as meal and oil from rapeseed, was 
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based on their economic value (Cederberg and Flysjoe 2004). 
Additionally, direct land use impacts were ascribed to the 
region in which the crop was most likely cultivated according 
to Cederberg and Flysjoe (2004) and for organic soy 
according to FiBL (2012). In this study, it was assumed that 
the land was occupied for one whole year for most crops, as in 
temperate latitudes only one fodder crop can be grown per 
year and oil palm fruit, meadows and pastures are cultivated 
permanently (Mila i Canals et al. 2013). Only after the harvest 
of soy another fodder crop can be grown in the same year. 
Thus, about 25 % double-cropping of soy was assumed (sug¬ 
gested in Dalgaard et al. 2008) which was in contrast to 
Cederberg and Flysjoe (2004) who did not account for any 
double-cropping. 


2.2 Land use assessment framework 

Following the framework of the United Nations Environment 
Programme/Society of Environmental Toxicology and Chem¬ 
istry Life Cycle Initiative (Mila i Canals et al. 2013; Koellner 
et al. 2013a, b), we distinguished two land use impacts: land 
occupation (using land) and land transformation (changing the 
land use) in this study (Mila i Canals et al. 2013; Koellner 
et al. 2013a, b). Permanent impacts (i.e. irreversible damages 
to ecosystems) were not considered here. Characterization 
factors of land occupation, CF 0cc? were calculated as the 
difference in relative plant species richness ( S re \) of a reference 
system and a land use type (LU) i per biome j (de Baan et al. 
2013; Koellner et al. 2013a, b). As globally available biodi¬ 
versity data focus on the species level and species numbers 
(Curran et al. 2011), we chose relative species richness of a 
certain land use type i compared to a reference as biodiversity 
indicator. N rel of the reference system (ref) was set to 100 %= 1 
(see Section 2.4.1). 


OF Occ^LUi,') — Srel,Ref)S re l,LUi,) ~ l _ ^re/,Zt/i,j (1) 

Natural vegetation was chosen as reference. In many Eu¬ 
ropean countries, only few anthropogenic undisturbed habitats 
remain, so they cannot be considered natural habitats (SOER 
Synthesis 2010). Consequently in the following, we will refer 
to the reference situations as "(sem^-natural", as proposed in 
de Baan et al. (2013). 

The numeric values of CF 0cc normally take on values 
between zero and one (expressing detrimental impacts), but 
can sometimes also be negative (expressing beneficial im¬ 
pacts, de Baan et al. 2013). The occupation impact can then 
be calculated as the product of CF 0cc? the cultivated area A 0cc 
and the duration of the occupation process t Qcc (de Baan et al. 
2013; Koellner et al. 2013a, b). 

Occupation impact LU = Aq cc x t( )cc x CFq C c (2) 


Transformation impacts are assessed accordingly with the 
transformed areaA Xrans (see Section 2.2) and the time a system 
requires to regain the ecosystem quality of the reference (4 eg ) 
after an anthropogenic disturbance (de Baan et al. 2013; 
Koellner et al. 2013a, b). 

Transformation impact Ref _> LU = A Trans x CF Trans 

A Trans ^ ^ ^ cc ^ ^ re S (^) 


The Biodiversity Damage Potential (BDP) of land use can 
then be calculated as the sum of the transformation impact and 
the occupation impact over all land use types i and biomes j 
(Koellner and Scholz 2007). 


BDP = ^Transformation Impact Re j_ 

A~ y Occupation Impact LU 


>LU, 


hj 



Table 1 General characteristics of the two farm types 



Organic 

Conventional 

Number of studied farms 

6 

9 

Arable land (ha farm *) 

63 

70 

Natural meadows (ha fann ' ) 

19 

11 

Number of cows per farm 

39 

65 

Stocking rate (LU a ha 1 year -1 ) 

0.9 

1.2 

Total purchased feed (kg cow -1 year -1 ) 

1,457 

2,951 

Milk yield (kg ECM b /cowxyear) 

9,400 

10,100 


a Lifestock unit: corresponds to one daily cow including a calf younger 
than 1 month, or six calves aged between 1 and 6 months, or three heifers 
older than half a year 

b Energy corrected milk 


2.3 Inventory analyses 

In the study of Cederberg and Flysjoe (2004), transformation 
impacts are not considered. Thus, we calculated the inventory 
data for transformed area as proposed by Mila i Canals et al. 
(2013). This approach only associates direct land transforma¬ 
tion with a fodder crop if (1) in its country of origin the 
harvested area of this specific crop increased in the last 
20 years and if additionally (2) the area of its land use type 
(i.e. arable land, permanent crops or meadows and pastures) 
increased. In case these two conditions applied, the 
transformed area for every occupied hectare and year was 
calculated by dividing the increase in land use type area over 
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the last 20 years by the current area of this land use 
type (as proposed in Mila i Canals et al. 2013). The 
type of land that was transformed was proportionally 
assigned based on the decrease in area of all land use 
types. Historical data on country-specific developments 
of agricultural production areas were extracted from the 
statistical database of the Food and Agriculture Organi¬ 
sation (FAOSTAT 2012a, b For grass pellets, data were 
extracted from Statistics Denmark (StatBank Denmark 
2012) as no information was given in the FAOSTAT 
database. For organic land use types, no specific infor¬ 
mation on land use change was available and standards 
of the International Federation of Organic Agriculture 
Movements only prohibit the conversion of areas with 
high conservation value 5 years before organic farming starts 
(IFOAM 2012). Therefore, we assumed organic farming to 
result in the same amount of transformed area per cultivated 
hectare as conventional farming. To account for yearly varia¬ 
tion in land use, 5-year averages around the years 1987 and 
2007 were calculated, respectively (as suggested in Mila i 
Canals et al. 2013). As data on forest area before 1990 were 
not available for all countries, pre-1990 areas were linearly 
extrapolated using data from the years 1990 to 2009. This 
seems a feasible estimate, as the forest conversion rate after 
1990 followed a linear relationship for the countries consid¬ 
ered in this study. Furthermore, within the timeframe in ques¬ 
tion (1985-2009) the same drivers were prevalent, as demon¬ 
strated in various studies (e.g. Feamside 2005; Gasparri et al. 
2008; McMorrow and Talip 2001). 

Forest is the only natural terrestrial land cover type 
included in the FAOSTAT database (2012a). Most Bra¬ 
zilian soy is cultivated on land where naturally savan¬ 
nahs prevail (Fearnside 2001). Savannahs still fall under 
the definition of '’forest” according to the FAOSTAT 
database, as they contain vegetation with trees exceeding 
5 m and a canopy cover denser than 10 %. Consequently, 
these data were seen as appropriate for estimating conversion 
from natural vegetation within the biome sub-/tropical grass-/ 
shrublands and savannahs. 

2.4 Impact assessment 

To account for spatial differences in ecological responses to 
land use, region-specific reference situations and regeneration 
times (see Eq. 3) were used and impacts were calculated 
separately for each biogeographic region (as suggested in de 
Baan et al. 2013). Biomes defined by Olson et al. (2001) and 
recommended to use for land use impacts in LCA on global 
scales by Koellner et al. (2013a, b) were used as biogeograph¬ 
ic units. CORINE Plus (elaborated by Koellner and Scholz 
2008) was used for land use classification as it allowed the 
discrimination of organic and conventional land use types 
(Koellner et al. 2013a, b). 


2.4.1 Biodiversity indicator 

The analysis of a biodiversity measure was restricted to vas¬ 
cular plants because data availability for organic land use 
types is relatively good for this taxon (Hole et al. 2005). 
Several data sources were combined here to quantify plant 
species richness of relevant land use types and reference 
situations: the Countryside Survey of the United Kingdom 
(Carey et al. 2007), the report "Biodiversity of Saxonian 
arable land" (Kreuter 2005), the biodiversity monitoring data 
of Switzerland (BDM 2004) and the databases used in 
Koellner (2000) and de Baan et al. (2013). As these data 
sources only contained little data on organic land uses or data 
from the biome sub-/tropical grass-/shrubland and savannahs, 
a literature search in the Web of Science database was 
performed. Overall, this search resulted in 66 studies (see 
Online Resource for a full bibliography) providing 111 data 
points for the different land use types and 53 data points for the 
reference situations in three different biomes of feedstock pro¬ 
duction for Swedish milk. As sampling area varied strongly 
among studies, sampled species richness (5) was standardized 
to an area (A) of 100 m using the transformed power model of 
the species-area relationship proposed in Kier et al. (2005): 

SlOOrn^ = S Samp ,e4 X (6) 

Sampled J 

where z is the species accumulation factor (see ESM, Table 1 
for results of each data point). For each biome, z values were 
obtained from Kier et al. (2005). 

2.4.2 Regeneration time 

The time required for the regeneration after human interven¬ 
tions stop is used for the calculation of transformation impacts 
(Eq. 3). To derive regeneration times of different land use 
types, a literature search of peer-reviewed studies was 
conducted. Only investigations located below 1000 m a.s.l. 
analysing regeneration of plant species richness after human 
disturbances (e.g. logging or agricultural use) were included. 
However, for the regeneration times needed for this study only 
very few data points (one to five per land use type and biome) 
could be found (Table 2). Therefore, alternative regeneration 
times based on Dobben et al. (1998) were calculated as well 
for a comparison. These regeneration times have been used in 
previous LCA studies on land use (e.g. Mila i Canals et al. 
2013), and are based on broad estimates of biomass regener¬ 
ation dependent on altitude and latitude after clear-cutting. To 
make this estimate comparable to the empirical data, biome- 
specific regeneration times were calculated depending on each 
biome’s area percentage (quantified with ArcGIS) in three 
different altitudinal and six different latitudinal zones. These 
values were considered to represent regeneration times of 
arable land (4 eg , arable,Dobben)* Thus, the regeneration time 
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Table 2 Regeneration times by land use type 



Based on Dobben 
et al. (1998) 

Leg (years) 

Mean of empirical 
studies 

Leg (years) 

n 

Source 

Biome sub-/tropical moist broadleaf forests 





Pasture/meadow 

45 

33 a 

5 

(Aide et al. 1995, 1996; Guariguata et al. 1997; 
Letcher and Chazdon 2009) 

Permanent crops 

42 

35 a 

3 

(Grau et al. 1997; Pascarella et al. 2000; 
Zimmermann et al. 2007) 

Biome temperate grass-/shrublands and savannahs 





Arable 

— 

45 b 

1 

(Scott and Morgan 2012) 

Biome sub-/tropical grass-/shmblands and savannahs 





Arable 

57 

— 




n number of data points per biome and land use type 
a All studies were conducted in Middle or South America 
b From Australia 

(4e g ,LUi,Dobben) of every land use type i was then calculated 
using the ratio of the CF 0cc per biome j (see Section 2.1) as 
follows: 

_ CF 0cCj Lu h j ^ . 

treg,LUi,j,Dobben ^ ^ reg,arablei,],Dobben {') 

C-L Occ ^arable 

It was assumed that the regeneration rate, defined as the 
change of CF 0cc , LUi with Leg? of each land use was constant 
within each biome. 

2.4.3 Statistical analysis 

No studies were found which had investigated the influence of 
organic and conventional farming practice and the species 
richness on comparable (semi-)natural references together. 
Therefore, we performed a resampling procedure with 1,000 
bootstrap samples (with replacement) to obtain several statis¬ 
tical parameters (e.g. median, quartiles, 95 % confidence 
intervals) of the relative species richness S TG \. For each land 
use type, bootstrapping combined all plant species richness 
data points randomly with (semi-) natural reference data 
points from the same biome. 

Analyses of variance (ANOVA) was applied to test differ¬ 
ences in means of S K \ among the factors land use type, biome and 
farming practice as well as their interactions. Prior to statistical 
analyses, bootstrapping results were log-transformed to improve 
fit to a normal distribution. 

As data did not satisfy the assumptions of the ANOVA (i.e. 
normal distribution and homogenous residual variance), ro¬ 
bustness of results was evaluated by additionally conducting 
the Kruskal-Wallis test to analyse the influence of the three 
factors (without interaction) on medians of A re] . The Mann- 
Whitney U test revealed significant differences in median S rel 
between pairs of land use types. All data were performed using 


R 2.13.0 (R Development Core Team 2011) and the R package 
“boot” (Angelo and Ripley 2012) for the bootstrapping. 

2.4.4 Biodiversity weighting factor 

Expressing species richness in relative instead of absolute 
terms allowed the comparison of land use impacts in various 
ecoregions with different data qualities (Kier et al. 2005; 
Koellner et al. 2013a, b). To account for differences in absolute 
species numbers and conservation value between ecoregions, a 
weighting system was applied, as proposed in Weidema and 
Lindeijer (2001) which was calculated separately for each 
ecoregion. Following global biodiversity prioritisation con¬ 
cepts, the biodiversity weighting factor (BWF) was based on 
absolute species richness , irreplaceability and vulnerability 
(see review by Brooks et al. 2006). 

The three indices to quantify biodiversity value of each 
ecoregion were calculated as follows. Absolute species rich¬ 
ness (S) was calculated as area-corrected total number of 
amphibian, reptile, mammal and bird species per ecoregion 
(derived from Kier et al. 2005). Irreplaceability was quantified 
as the area-corrected number of strict endemic species of 
amphibians, reptiles, mammals and birds (EndS). For ende¬ 
mism, these are the only taxonomic groups where data per 
ecoregion are available. For consistency, the same selection of 
taxonomic groups was also chosen for species richness, and 
data on plants were excluded. Vulnerability was expressed as 
the 'Conservation Risk Index' (CRI), which is calculated as the 
ratio of converted ecoregion area in per cent to protected 
ecoregion area in per cent (developed by Hoekstra et al. 
2005). The latter concept assumes that the more area is occu¬ 
pied the more damaging an occupation or transformation will 
be for the remaining ecosystem (Koellner 2000). To prevent 
division by zero, all values below 1 % were set to 1 % as such 
low habitat proportions are likely below the threshold necessary 
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for long-term species conservation (Swift and Hannon 2010). 
The same procedure was performed for converted area as the 
impact of converting 1 % of a habitat is assumed to be negli¬ 
gible in terms of species conservation. 

Data to calculate all indices were extracted from the Terres¬ 
trial Ecoregions Base Global Dataset (Olson et al. 2001). As 
suggested by Weidema and Lindeijer (2001), these three mea¬ 
sures of conservation prioritisation were linked by multiplica¬ 
tion for each ecoregion k. To give each measure equal weight in 
the multiplication, they were all normalised (expressed with an 
asterix * in equation) to range between 1 and 10 (Table 3): 


BWFv = 


'4x9 
S 


(EndS k x 9 


max 


+ n x "' +1 


EndS 


CRI k x 9 

x I —^-1- 1 


CRI 


max 


= S k * x Ends** X CRI k 


( 8 ) 


Overall, the weighted BDP was calculated as follows: 


weighted BDP k 


(E TI Ref^LU,, k + E OIlU <*) 


X BWF k 



3 Results 

3.1 Characterization factors for biodiversity damage potential 

For almost all land use types, medians of CF 0cc were signif¬ 
icantly greater than zero, indicating detrimental impacts on 


Table 3 Results of normalisations of species richness (S), endemic 
species richness (EndS) and Conservation Risk Index (CRI) and their 
product, the BWF per ecoregion 


Biome 

Ecoregion 

S a 

EndS a 

CRI a 

BWF 

SubVtropical moist 
broadleaf forest 

Peninsular 

Malaysian rain 
forests 

4.6 

1.5 

3.0 

18.1 

Temperate broadleaf & 
mixed forests 

Atlantic mixed 
forests 

2.4 

1.0 

7.8 

18.6 


Baltic mixed 
forests 

2.4 

1.0 

5.4 

12.8 


Sarmatic mixed 
forests 

2.1 

1.0 

1.6 

3.4 

Sub-/tropical grass-/ 
shrublands and 
savannahs 

Cerrado 

3.3 

1.8 

5.4 

31.7 


Sources: Olson et al. (2001) and own calculations. A full list of BWF per 
ecoregion can be found in Online resource Table 6 

a Results of normalisations 


biodiversity (see Figs. 1 and 2, and ESM Fig. 1, Table 2 and 
Table 6). Only the median of organic pasture/meadow was not 
significantly different from zero and showed a trend towards a 
negative value (i.e. beneficial impact). Overall, meadows/ 
pastures were least harmful to biodiversity, followed by per¬ 
manent crops. Arable land had the highest impact on biodiver¬ 
sity, which was most pronounced in the biome sub-/tropical 
moist broadleaf forests. A significant effect of land use type, 
biome and farming practice on CF 0cc was shown with ANOVA 
and Kruskal-Wallis test (results not shown). The ANOVA also 
revealed that the influence of land use type on CF 0cc was 
significantly affected by biome and farming practice. 

In all studied land use types and biomes, organic farming 
practices were less detrimental for biodiversity than conven¬ 
tional farming practices. Medians of CF 0cc of conventional 
compared to organic pasture/meadow showed a 30 % higher 
potentially disappeared fraction (PDF) of species per square 
meter and year. In the median, 45 % less species were found 
on conventional arable land in the biome temperate broadleaf 
and mixed forests than on the organic equivalent. In the biome 
sub-/tropical grass-/shrublands and savannahs, this difference 
amounted to about 40 %. 

3.2 Impact of occupation 

Although organic milk required about double the amount of 
agricultural land to produce 1 kg of milk (Fig. 3 a), the occu¬ 
pation impact of organic milk was only half the one of con¬ 
ventional milk (Fig. 3b). CF 0cc of organic land use types were 
always considerably lower than the conventional ones thus 
leading to smaller occupation impacts. In addition, the differ¬ 
ent composition of the feedstock, with larger shares of rough- 
age feed and grazing for the organic cows and larger shares of 
concentrate feed for the conventional cows, considerably 
influenced the result (ESM, Table 3). 

Arable land and more precisely grains and legumes/ 
temporary grasses in the biome temperate broadleaf and 
mixed forest largely dominated the impact for both farming 
practices (Fig. 4). For the fodder crop category legumes/ 
temporary grasses 100 and 80 % of the impact was due to 
ley for conventional and organic milk, respectively. Both 
conventional and organic farms cultivated ley on about 60 % 
of their on-farm arable land (Cederberg and Flysjoe 2004). 
Grains accounted for more than 75 % in the organic and about 
50 % in the conventional concentrate feed (ESM, Table 3). 
This explains the large land use and occupation impact allo¬ 
cated to this crop type. In contrast, rapeseed, sugar beet and 
palm oil have relatively small shares in the occupation impact 
of conventional milk, although they constitute about 37 % of 
the concentrate feed. Due to their form as co-products of the 
sugar/oil/starch production, only a small part of the environ¬ 
mental burden was allocated to them, minimizing their impact. 
Although organic milk required a relatively large area of 
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Fig. 1 Box and whisker plot of bootstrapped characterization factors of 
occupation (CF 0cc ) for each land use type, farming practice and biome. 
TempBMF temperate broadleaf and mixed forests; TropGL sub-/tropical 
grass-/shrublands and savannahs; TropMBF sub-Aropical moist broadleaf 
forests. Asterisks mark significant differences in medians compared to the 
reference which is defined as CF Occ = 0 (Mann-Whitney U test,/?<0.05). 


n is the total number of data points used for the bootstrap of CF 0cc ; the 
first number stands for the number of species richness data points per land 
use type and the second for the data points per reference. 77 = 14 x 33, e.g. 
means that 14 data points on the species richness of the land use type were 
combined with 33 data points for the reference situation, which resulted in 
a total of 462 data points for CF 0cc 


pastures and meadows, the occupation impact was slightly 
negative (i.e. beneficial) due to their negative median CF 0cc 
value. In the organic feedstock, some conventionally pro¬ 
duced grains are included. This is the reason why grains and 
legumes have a similar occupation impact (Fig. 3b), although 
larger areas are required for legume production than for grain 
production (Fig. 3a). 


Biome TropGL 



Transformation from: Forest to arable Permanent crops Forest to 

to arable permanent crops 

n=15x6x1 n=2x6x1 n=12x14x3 

Fig. 2 Box and whisker plot of characterization factors of transformation 
(CF Trans ) for each farming practice, biome and land use type regenerating 
to after human abandonment. TropGL sub-/tropical grass-/shmblands and 
savannahs; TropMBF sub-/tropical moist broadleaf forests. Shown is only 
the variation due to differences in species richness between data points. 
See Online Resource Fig. 1 for variation in species richness per land use 
type and biome 


3.3 Impact of transformation 

Conventional milk had a more than threefold larger transfor¬ 
mation impact than the same amount of organic milk, regard¬ 
less of the method used for estimating regeneration times 
(Fig. 5). In this study, soy and palm oil were the only fodder 
crops originating from countries, where in the last 20 years the 
harvested crop area as well as the area of the corresponding 
land use type (i.e. arable land, pastures/meadows or perma¬ 
nent crops) increased. Soy was responsible for the majority of 
the transformation impact (99 % for conventional milk and 
98 % for organic milk). In contrast, palm oil caused only a 
very small impact because palm kernel expels have only little 
economic value and the quantities in the concentrate feeds 
were small (3 % for conventional milk and 4 % for organic 
milk in concentrates). The land use change from natural 
vegetation to each harvested hectare of crop was the same 
for the conventional soy from Brazil (0.20 ha converted/ 
harvested ha per year) and palm oil from Malaysia (0.20 ha 
converted/harvested ha per year), but slightly smaller for the 
organic soy from Argentina (0.14 ha converted/harvested ha 
per year). Argentina was the only country in the milk value 
chain, where one agricultural land use type, permanent crops, 
was converted into another agricultural land use type (ESM, 
Table 4). Thus the transformation impact of Argentinian soy is 
lower than Brazilian soy as the regeneration time of arable 
land to permanent crops is smaller than from arable land to 
secondary vegetation. Also, the much higher yields of con¬ 
ventional compared to organic soy were counterbalanced by 
the lower CF Trans of organic arable land in this biome. As a 
consequence, the difference in the transformation impact 
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Fig. 3 Contribution to a land 
occupation and b occupation 
impact of different fodder crop 
types per biome TempBMF 
temperate broadleaf and mixed 
forests, TropGL sub-Aropical 
grass-/shrublands and savannahs 
and TropMBF sub-/tropical moist 
broadleaf forests, calculated with 
non-weighted characterization 
factors for biodiversity impact, 
which is expressed as potentially 
PDF of species per FU 



<D 

CL 

CO 

CM 

E 




between organic and conventional milk can largely be 
explained by the fact that conventional cows are fed about 
three times more soy products than organic cows to produce 
the same amount of milk. 

3.4 Total biodiversity damage potential 

The total BDP is about 60 % lower for organic compared to 
conventional milk. The share of the transformation impact on 
the BDP amounts to about 30 and 20 % for conventional and 
organic milk, respectively (Fig. 6a). From all feedstock, soy 
showed the largest impacts on biodiversity. Although soymeal 
comprised only about 10 % of the conventional concentrate 


feed, it caused almost 40 % of the BDP. In the case of organic 
milk, soy constituted only 4 % in the organic concentrate feed, 
but was responsible for 25 % of the BDP. 

To account for differences in absolute species numbers and 
conservation value between ecoregions, we calculated a biodi¬ 
versity weighing factor. For the ecoregions which are relevant 
for the production of Swedish milk, the normalized weighing 
factor ranged between 4.3 for the Sarmatic mixed forest and 
31.7 for the Cerrado (Table 3). The ratio in the weighted BDP 
of organic to conventional milk was slightly decreased from 43 
to 36 %. The impact of rapeseed from the Baltic mixed forest, 
which was negligible without the weighting system, became 
obvious. After weighting, soy largely dominated results of the 
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Biome TempBMF BiomeTropGL Biome TropMBF 



Pasture/ Arable Arable Permanent Total 
meadow crops 


Fig. 4 Occupation impact of conventional and organic milk illustrated 
per land use type and biome TempBMF temperate broadleaf and mixed 
forests, TropGL sub-/tropical grass-/shrublands and savannahs and 
TropMBF sub-/tropical moist broadleaf forests, calculated with non- 
weighted characterization factors for biodiversity impact, expressed as 
PDF of species per FU 

BDP. This is due to the relatively high weighting factor of the 
Cerrado ecoregion where most of the soy in this study was 
cultivated. Fifty percent of the Cerrado has been converted, 
whereas only 1.1 % of its area is protected, resulting in a high 
CRI value. 



Fig. 5 Transformation impact calculated with regeneration times based 
on a empirical data and b expert estimates of biomass regeneration by 
Dobben et al. (1998), calculated with non-weighted characterization 
factors for biodiversity impact, expressed as PDF of species per FU. 
TropGL sub-/tropical grass-/shrublands and savannahs; TropMBF sub-/ 
tropical moist broadleaf forests 


4 Discussion 

For some, organic agriculture is seen as solution for environ¬ 
mentally friendly food production, while for others, it is an 
inefficient and resource-intensive production system. To bet¬ 
ter understand the involved trade-offs, quantitative compari¬ 
sons are needed. To our knowledge, this is the first study 
trying to quantify and compare the land use impacts on 
biodiversity of organic and conventional farming on a product 
level, including impacts of imported feedstock. While organic 
milk required about double the area than conventional milk, 
the impacts on biodiversity were less than half for organic 
compared to conventional milk. The impacts on biodiversity 
per organically farmed area were considerably smaller for 
conventional farming and less fodder crops originated from 
sub-/tropical countries where forest conversion to agricultural 
land uses is prevalent. These reduced impacts of organic milk 
could by far outweigh the larger area requirements. 

Trying to express such complex issues as global land use 
impacts on biodiversity in one number is inherently accom¬ 
panied by large natural variation and uncertainties. In the 
following, we want to discuss the sensitivity and robustness 
of results. In addition, we will discuss methodological choices 
made in this study, such as how to quantify biodiversity 
impact and the amount of transformed area, the selection of 
the indicator for regeneration times and the appropriateness of 
the weighing factor. Finally, we will reflect on the transfer- 
ability of results of this study to milk produced in other 
countries. 

4.1 Sensitivity analyses 

This paper combined data from a Swedish LCA study, global 
literature on plant species richness and FAO statistics. All 
databases exhibit diverse sources of uncertainty. Conducting 
a full Monte Carlo uncertainty analyses was not possible in 
this study as information on uncertainty or variability was not 
available for all parameters. Therefore, to assess how the 
uncertainty of selected parameters influenced the results, a 
sensitivity analysis was performed (see Fig. 7 for an 
overview). 

Inter-farm variation is rarely quantified in LCA studies on 
agricultural products, but can be considerable. Land use in¬ 
ventory data for this case study came from a survey of nine 
conventional and six organic dairy farms from the same region 
within Sweden. For these farms, the variability within the land 
occupation of both farming practices is remarkable; the farm 
with highest land use needs twice the area per functional unit 
than the one with lowest, in both farming practices, respec¬ 
tively. To test the robustness of the results, we compared milk 
from the organic farm with the largest land use and milk from 
the conventional farm with the smallest land use (see Fig. 7 
scenario SI). In this latter case, the organic farm had a five 
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Fig. 6 Biodiversity damage 
potential (sum of occupation and 
transformation impact) calculated 
a without and b with the 
biodiversity weighting factor per 
ecoregion, expressed as PDF of 
species per FU. Be aware of the 
different scales of the axes 
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times larger land use than the conventional farm. This repre¬ 
sents an extreme case for Europe, as in the published studies 
of other countries differences in land requirements between 
organic and conventional milk were less pronounced 
(Thomassen et al. 2008). In this case, the direct land occupa¬ 
tion impact of organic milk from the farm with the largest land 
use was slightly higher than conventional milk from the farm 
with the smallest land use. Understanding the factors causing 
this inter-farm variability would be an important topic for 
future research as it might offer optimisation potential for both 
farming practices. 

In agricultural LCA studies, yields are often one of the 
determining factors. In this current study, organic compared to 



Fig. 7 Sensitivity of biodiversity damage potential (sum of occupation 
and transformation impact) per FU for scenario SO to S5 {SO standard 
scenario used in this study; SI comparing organic farm with highest land 
use with conventional farm with lowest land use; S2 two times lower 
organic yields than used here; S3 South American soybean meal substitut¬ 
ed with North American soybean meal; S4 ley classified as pastures and 
meadows instead of arable land; S5 using third quartiles instead of 
medians of characterization factors of occupation) 


conventional yields were estimated to be 20 to 50 % lower per 
crop. Compared to yield differences ranging between 5 and 
34 % found in a global meta-study (Seufert et al. 2012), our 
yield estimates are conservative. Even if the yields of organic 
crops were on average about two times lower than adopted in 
this study (assuming a constant composition and origin of 
feedstock for both alternatives), the organic milk would still 
be favourable in terms of impacts on biodiversity (see Fig. 7 
scenario S2). Purchased feedstock is often sourced from var¬ 
ious countries. As we did not have data for all world regions 
on plant species richness of organic and conventional agricul¬ 
ture, we could not directly test how the results would change 
assuming a different origin. As soy meal was responsible for 
about 25 % of the total BDP of organic and about 40 % of 
conventional milk, the total direct land use impacts could 
change significantly, when sourcing soy meal from another 
country, such as the USA, China or India. According to the 
applied methodology, in those countries, no transformation 
impact would have been allocated to soy, as in these countries 
the area of arable land remained stable for the last two decades 
(FAOSTAT 2012a) therefore decreasing the BDP (see Fig. 7 
scenario 3). 

Uncertainty also lies in the classification of land uses. Feys 
cultivated with temporary grasses and legumes grown on 
areas set aside were classified as arable land in this study. 
The reason for this decision was that leys were cultivated on 
average for 3 to 8 years in a crop rotation system with regular 
ploughing. According to Koellner et al. (2013a, b), the latter is 
an important differentiating factor from pastures/meadows in 
the land use and land cover classification used in this study. A 
Norwegian study found that the species richness of leys 
reached levels comparable to semi-natural grassland not be¬ 
fore 30 years after last ploughing (Austrheim and Olsson 
1999). However, ley grasslands might have a higher species 
richness than average arable fields cultivated as monocultures 
(Andreasen and Stryhn 2008). Thus, we performed a 
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sensitivity analysis classifying all the ley area as meadow/ 
pasture. This decreased the occupation impact of conventional 
and organic milk about 22 and 46 %, respectively, and resulted 
in even more pronounced differences between organic and 
conventional milk (see Fig. 7 scenario S4). However, up to 
now almost no studies were conducted on the plant species 
richness of leys. Thus, more research is needed to improve 
reliability of biodiversity impact assessments dealing with 
land use on leys. 

An in-depth uncertainty assessment of direct land use 
impacts on biodiversity of globally traded cocoa products 
showed that the CF 0cc were by far the largest contributor to 
variance (Mutel et al. 2013). In this study, the variance of 
CFqcc was also considerable, even ranging from positive to 
negative impacts for pasture and arable land in the temperate 
biome and for permanent crops in the tropical biome (see 
Fig. 1). Part of this variation might be due to the fact that 
different data sources had to be used for the plant species 
richness of the different land use types and of the reference. 
Studies directly comparing the species richness of organic 
agriculture and a reference within the same region were not 
available. Consequently, data points of land use types had to 
be randomly combined with (semi-) natural reference data 
points from the same biome to generate CF 0cc > possibly 
resulting in plant species richness data from species-poor 
regions being combined with (semi-) natural references from 
species-rich regions in the same biome. 

We tested the sensitivity to CF 0cc by recalculating the BDP 
results based on the third quartiles of CF 0cc Using the third 
quartiles instead of medians of CF 0cc increased the BDP for 
conventional and organic milk by about 120 and 230 %, 
respectively (see Fig. 7 Scenario S5). In this case, organic 
milk still had only half of the impact of conventional milk. 
When using the first quartile of CF 0cc for conventional land 
use types and the third quartile of CF 0cc for organic land use 
types, the result changes and the BDP of conventional milk 
would be lower than the one of organic milk. Thus, future 
research should aim to decrease the uncertainty of CF 0cc by 
possibly combining data points on species richness of land use 
types with data points on (semi-) natural references from the 
same region. However, a large variation will possibly still 
remain, as biodiversity shows a strong natural heterogeneity. 

4.2 Quantification of biodiversity impacts 

The strong difference in plant species richness between or¬ 
ganic and conventional farmed land found in this study has 
also been found in other studies (e.g. Bengtsson et al. 2005; 
Fuller et al. 2005; Gomiero et al. 2011). A meta-analysis by 
Hole et al. (2005) on the influence of organic agriculture on 
biodiversity reported that differences between the two farming 
practices were especially pronounced for plants, birds and 
insects. However, for other taxonomic groups, such as 


earthworms or beetles, the influence of organic farming on 
species richness was less clear, as studies reported both pos¬ 
itive and negative effects (Hole et al. 2005). Another meta¬ 
analysis by Bengtsson et al. (2005) found that the effect of 
organic agriculture was largest on studies on the plot scale, 
whereas on the farm scale results were still significant, how¬ 
ever less pronounced. This might suggest that the kind of data 
that were chosen in this study (i.e. plant species richness on 
field scale or smaller) made results more clear than it could 
have been in the case if all organism groups had been studied. 
However, plants were the only taxon for which sufficient data 
were available for all relevant biomes and land use types. 
Moreover, two studies (Haas et al. 2001; Schader et al. 
2010), which compared biodiversity impacts of organic and 
conventional milk production systems on farm scale, also 
found lower impacts of organic dairy farms. However, their 
approach did not allow the quantification of impacts related to 
off-farm land use which this study showed to be responsible 
for a considerable part of the BDP of milk. 

Another factor which seems to have an influence on the 
species richness of agricultural land is the surrounding land¬ 
scape. In their meta-analyses, Bengtsson et al. (2005) con¬ 
cluded that positive effects of organic farming will more likely 
occur in intensively used agricultural land, but not inevitably 
in small-scaled mosaic landscapes with a high share of natural 
vegetation. Thus, the results found here seem to be more 
representative for the highly agriculturally used landscapes 
found in middle Europe. 

4.3 Calculation of transformed area 

The quantification of area transformed for land use is largely 
based on FAO statistics. The database has the advantage of 
being easily applicable and consistent across countries (Mila i 
Canals et al. 2013). However, one has to be aware of several 
shortcomings of the database in this context. First, the lack of 
forest area data for the last 20 years poses a problem for 
studies like this one. Calculations can either be based on 
extrapolations of existing data or periods shorter than 20 years. 
This influences the results (compare with Mila i Canals et al. 
2013). A second shortcoming is that the only natural terrestrial 
vegetation type included is '’forest'’ (see Section 2.2 for the 
definition). All areas which are not classified as agricultural 
land or forest area are compiled in "other land" which e.g. 
sums up to 37 % of the country area in the case of Argentina. 
Hence, the assessment of conversion from more grassland- 
type vegetation is difficult. This study tried to overcome this 
problem by assuming that the entire increase of the land use 
type of a crop was due to direct land transformation from 
natural vegetation in case no other agricultural land use type 
decreased in area during the respective time period in the 
country. Furthermore, uncertainty of the FAO database is high 
as many data are "FAO estimates" or "unofficial figures". 
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Thus, in some cases it could be appropriate to compare FAO 
data with national inventory data, as performed for land use 
estimates by Ramankutty et al. (2008). Moreover, the method 
does not consider indirect land use change (iLUC) in other 
countries (Schmidt 2008b; Mila i Canals et al. 2013) which 
might be relevant in comparing different farming practices. It 
is possible that organic products have a higher iLUC as their 
yields are lower than conventional products (Seufert et al. 
2012). However, in the scope of this study it was not possible 
to model worldwide socio-economic impacts of lower organic 
yields compared to conventional fodder crops on iLUC in 
other countries. Including iLUC into environmental impact 
assessments has only developed in the last years and was 
mainly limited to calculations of carbon emissions due to 
biofuel cultivation (e.g. Kloverpris and Mueller 2013). Meth¬ 
odological considerations and assumptions concerning iLUC 
are highly debated (Kloverpris and Mueller 2013; Hertel et al. 
2010; Mathews and Tan 2009) and estimates vary consider¬ 
ably (Plevin et al. 2010). Nevertheless, it is expected that 
parameter values used to model iLUC of biofuels are not 
comparable to the ones that would be needed to model iLUC 
of organic crops. One parameter which is different is price: 
organic milk is about 15-83 % more expensive in the Euro¬ 
pean Union (Austrian Federal Ministry of Agriculture, For¬ 
estry, Environment and Water Management 2011) which is 
known to decrease milk consumption (Audreyeva et al. 2010) 
and thus expected iLUC. However, the socio-economic 
cause-effect chain of converting agricultural land to organic 
land is complex and potential impacts on land management in 
other countries are barely understood. Moreover, not enough 
data to calculate CF 0cc or regeneration times for other biomes 
where natural ecosystems are converted for agricultural use 
were available. Nevertheless, including iLUC into the assess¬ 
ment of greenhouse gas emissions of biofuels was shown to 
cancel out their benefit (e.g. Hertel et al. 2010). Therefore, in 
future comparisons of different farming practices related to 
land use change should include this aspect. 

4.4 Indicator for regeneration time 

Up to now, biomass estimates (e.g. Mila i Canals et al. 2013; 
Schmidt 2008a) were used in LCIA as indicator for the regen¬ 
eration of biodiversity as a whole, in contrast this study also 
used empirical data of plant species richness. The similarity of 
the results for transformation impacts (see Fig. 5) calculated 
with empirical data in comparison with biomass estimates by 
Dobben et al. (1998) is astonishing. Dobben et al. (1998)'s 
estimates depend only on latitude and altitude, whereas the 
regeneration of species richness is also modified by factors 
such as the natural disturbance regime (Brown and Lugo 
1990), soil nutrients or the proximity to patches of natural 
vegetation in the surrounding landscape (Prach and 
Rehounkova 2006). But what takes much longer than the 


regeneration of biomass or vascular plant species richness 
are the regeneration of species composition (Aide et al. 
2000; Pascarella et al. 2000; Dunn 2004), species richness of 
other taxa (Barlow et al. 2007; Berry et al. 2010) or endemic 
species richness (Liebsch et al. 2008). A recent meta-analysis 
of empirical studies showed that primary tropical forests are 
irreplaceable in terms of biodiversity (Gibson et al. 2011). A 
practical solution for this problem in LCA applications could 
be to assume a very long regeneration phase (Koellner and 
Scholz 2007). Koellner et al. (2013a, b) recommended to limit 
the modelling period in land use impact assessments to 
500 years. Incorporating such long regeneration times would 
increase the share of transformation impacts on the overall 
BDP considerably. Ultimately the choice of indicator depends 
on the aim of biodiversity LCIA (see Michelsen 2011). 
Depending on whether we focus on ecosystem services, func¬ 
tional diversity should be protected (see also de Baan et al. 
2013) which correlates well with plant species richness 
(Petchey et al. 2004). However, if we focus on the conserva¬ 
tion or the potential prospective economic value (e.g. medic¬ 
inal plants) of biodiversity, endemic species richness or spe¬ 
cies composition of primary habitat should be chosen. A meta¬ 
analysis on empirical regeneration times of species richness 
for different world regions and taxa will become available 
soon (Curran et al. in press) which will improve reliability of 
results as uncertainty of regeneration times used in this study 
is high due to few available data points (see Table 2). Regen¬ 
eration times in this study (see Table 2) of tropical biomes are 
not directly comparable to the numbers found for temperate 
biomes as the compared land use types for the biomes differ. 
Regeneration of species richness on arable land is expected to 
take much longer than on pastures/meadows or permanent 
crops as the human intervention is much more intense due to 
regular ploughing. 

4.5 Biodiversity weighting system 

To account for differences in absolute species numbers as well 
as conservation value between ecoregions, a weighting factor 
was introduced, as recommended in de Baan et al. (2013). The 
factor proposed here is easy to apply and consistent over 
ecoregions on a global level. However, up to now only infor¬ 
mation for amphibians, reptiles, mammals and birds could be 
included. Thus, as soon as data on endemic plant species or on 
other important species groups such as arthropods become 
available, additional taxa should be included in the calculation 
of a weighting system to improve the information value. 

Here, vulnerability was assessed as the ratio of converted to 
protected ecoregion area. Yet, converted areas can also be a 
main target of nature conservation for example the species- 
rich grasslands developed after deforestation and subsequent 
extensive grazing and mowing in Europe (Poschlod and Wal¬ 
lis de Vries 2002). Therefore, the vulnerability assessment 
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might be improved by calculating the ratio of natural and 
extensive land uses to protected ecoregion area. Including 
the areas of extensively used land might have resulted in a 
higher biodiversity weighting factor of the Sarmatic mixed 
forest. Currently, data on areas of extensive or traditional land 
uses are not available on ecoregion scales. Nevertheless, re¬ 
sults of the introduced weighting scheme in this study showed 
that the Cerrado exhibits a considerably more valuable biodi¬ 
versity than European ecoregions. Accounting for this remark¬ 
able richness is essential in the assessment of biodiversity 
impacts in LCAs. 

4.6 Transferability to other temperate countries 

Results for occupation impacts of milk found in this study 
cannot be transferred easily to milk produced in other coun¬ 
tries. When comparing the land use inventories of various 
European milk LCA studies, the variation in allocation pro¬ 
cedures and assumptions is posing a difficulty (De Boer 2003; 
De Vries and De Boer 2010; Yan et al. 2011). The sources of 
feedstock vary, where some studied farms source concentrate 
feed components from countries such as India (De Boer et al. 
2013) or Australia (Thomassen et al. 2008) The biggest dif¬ 
ference across countries seems to be their share of pastures/ 
meadows on the total land use. As CF 0c c of pasture/meadow 
are relatively low, high proportions of this land use decrease 
the occupation impact. In Sweden, land use relies to only 12 % 
for conventional and 20 % for organic milk on permanent 
pastures/meadows (calculation based on data in Cederberg 
and Flysjoe 2004). These percentages do not include the ley 
areas which seem to be special for Sweden but are classified as 
arable due to their regular ploughing. Thus, land use for 
Swedish milk production is based less heavily on meadows/ 
pastures compared to Belgium with 42 % for zero-grazing (i.e. 
no grazing on pastures but roughage from meadows fed) 
conventional milk (Meul et al. 2012), Germany (Mueller- 
Lindenlauf et al. 2010) with 45-89 % for organic milk or the 
Netherlands (Thomassen et al. 2008) with 43 % for conven¬ 
tional and 54 % for organic milk. The highest share of per¬ 
manent pastures/meadows on the overall land use for milk 
seems to be found in New Zealand (Basset-Mens et al. 2009) 
where 95 % of the land use for high-input milk and 100 % of 
the land use for low-input milk (i.e. no purchased feed outside 
the farm) is based on permanent mixed pastures. However, it 
has to be investigated whether CF 0c c generated in this study 
are transferable to other European countries or regions in the 
biome temperate broadleaf and mixed forest (e.g. New 
Zealand, China or USA) or to which extent other factors such 
as agricultural history or sensitivity of native species to agri¬ 
culture cause significantly different results. 

Transformation impacts found in this study seem to be 
transferable to milk produced in other European countries. 
Most middle and northern European dairy farms source their 


feed largely from countries where the area of agricultural 
production remained stable for decades (FAOSTAT 2012a). 
Therefore, the share of soy meal (from South America) and 
palm kernel expels (from South-East Asia) in the concentrate 
dairy feed might be a good proxy for transformation impacts. 
The percentage and origin of soy and palm kernel expels in 
conventional concentrate feed in this case study (about 10 % 
of 0.32 kg concentrate feed per FU) are comparable to other 
European countries. A Dutch study listed 12 % Brazilian soy 
of 0.26 kg concentrate feed per FU (de Boer et al. 2013), a 
French study 13.8 % soy from Brazil (Lehuger et al. 2009) and 
in the regional estimate of the FAO (2010) 13 %. In contrast, 
percentages in organic concentrate feed (5 %) in this study 
seem to represent values of the upper range. In Austria 
(Hoertenhuber et al. 2010), no soy meal or palm kernel expels 
were included in organic concentrates and almost none (ac¬ 
counting for economic allocation) in the Netherlands 
(Thomassen et al. 2008). Thus, transformation impact results 
calculated for conventional milk seem to be representative for 
other European countries, whereas results for organic milk 
seem to demonstrate values on the upper margin. Finally, the 
statement that organic milk has lower transformation impacts 
than conventional milk seems to be valid for the studies 
mentioned here. 


5 Conclusions 

Although organic land requires about double the area than 
conventional milk, the direct impacts on biodiversity were less 
than half. This illustrates the importance of differentiating 
CFqcc depending on the land use intensity (e.g. organic versus 
conventional). However, as including indirect land use change 
in assessments of greenhouse gas emissions of biofuels was 
shown to be crucial, future research on the biodiversity impact 
of organic products should study socio-economic effects of 
lower organic compared to conventional yields on indirect 
land use change. For an application to other studies, a more 
complete data set on species richness of organic and conven¬ 
tional farms needs to be compiled, to cover more world 
regions. However, for studies concerned with similar coun¬ 
tries of origin, the presented CFs can be used to approximate 
impacts, considering the underlying uncertainties as discussed 
above. To increase validity of results, more research is re¬ 
quired especially on the impact of land use on biodiversity 
compared to a (semi-) natural reference in grassland biomes, 
which are currently gaining in agricultural importance. Future 
research should assess the relevance of a more fine-scaled 
regionalization and develop concepts on incorporating aspects 
such as type of surrounding landscape or history of land use, 
which have a considerable influence in some regions 
(Bengtsson et al. 2005). 
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In the present study, we just compared the direct land use 
impacts to the endpoint Biodiversity Damage Potential 
(Koellner et al. 2013a, b) and did not consider other 
midpoint-oriented impact indicators (e.g. climate regulation, 
biotic production, freshwater regulation). To better understand 
the role of organic agriculture in an environmentally friendly 
solution to feed the world, other environmental as well as 
socio-economic impacts should be considered. As was shown 
in this study, further methodological developments might be 
needed to capture the differences between organic and con¬ 
ventional farming systems within LCA. 

Another finding of this study is that meadows and 
pastures are the most biodiversity-friendly feedstock for 
milk production. Meadows and pastures generally will 
have a higher species richness than arable land, assuming 
stocking rates are not detrimental to biodiversity. The 
larger the share of pastures and meadows on the land 
occupation of milk, the lower the impact to biodiversity 
will be. As a consequence, politics should create incen¬ 
tives for farmers to maintain pastures and meadows, if 
they contain high species richness. Further improvement 
potential exists by sourcing soy from countries without 
deforestation or to replace soy with another, more 
biodiversity-friendly crop. Results found here also stress 
the importance of subsidies for organic agriculture as this 
type of farming practice makes an important contribution 
to the maintenance of species richness in the agricultural 
landscape. 
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